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Abstract 

The nucleation of a lamellar phase from a supercooled homogeneous phase in a fluc- 
tuation driven first-order transition is studied, based on a phenomenological free energy 
due to Brazovskii. The absence of phase coexistence in the corresponding mean-field 
approximation makes application of the standard droplet theory of homogeneous nucle- 
ation problematic. A self-consistent coarse-graining procedure is introduced to overcome 
this difficulty, and the barrier height for nucleation of a critical droplet is estimated in 
the weak-coupling limit. Contrary to earlier estimates the critical droplet shape is 
shown to be anisotropic in general. Some effects of distortions and defects in the lamel- 
lar structure are considered and are shown to affect the critical droplet only very near 
coexistence where the probability of nucleation vanishes. The coarse-graining procedure 
introduced here follows from a novel application of the momentum-shell renormalization 
group method to systems in the Brazovskii class. Possible applications of the theory 
to the microphase separation transition in diblock copolymers and to Rayleigh-Benard 
convection are briefly discussed. 
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I Introduction 



In 1975 Brazovskii [1] investigated isotropic or nearly isotropic systems in which the fluctuation 
spectrum had a minimum at a nonzero wavevector |q| = q = q , represented by a hypersphere in 
ti-dimensional reciprocal space. He found that in the symmetric case where mean-field theory would 
predict a continuous transition to a periodic ordered state with spatial period 2ir/qo, fluctuation ef- 
fects lead to a discontinuous or first-order transition. In the case of weak fluctuations (small noise), 
the transition point is close to its mean-field value and the self-consistent Hartree approximation 
employed by Brazovskii could be justified as the leading term in a systematic perturbation expan- 
sion. The physical origin of the effect, which is essentially independent of system dimensionality 
for d > 2, lies in the large phase space for one-dimensional fluctuations in the direction transverse 
to the hypersphere. 

Although initially proposed to describe weakly anisotropic antiferromagnets and cholesteric 
liquid crystals, the Brazovskii model was subsequently shown to apply to the nematic to smectic-C 
transition in liquid crystals [2] , to pion condensates in neutron stars [3] , to the onset of Rayleigh- 
Benard convection [4,5] and most notably to the microphase separation transition in symmetric 
diblock copolymers [6-9]. It is in this last system that experimental confirmation of the theory 
was achieved [8] by quantitatively estimating the mean-field parameters and showing that the 
observed first-order transition was inconsistent with the predictions of mean-field theory, but was 
well described by the Brazovskii theory. 

Given this measure of success it is natural to ask about the lifetime of the metastable isotropic 
phase as the system is cooled below the thermodynamic transition point [10]. From a theoretical 
point of view the problem is formulated by classical homogeneous nucleation theory [11] as well 
as by more sophisticated field-theoretic approaches [11, 12], in terms of the free energy barrier 
for creation of a droplet of ordered phase inside the disordered phase. The critical droplet is the 
one which balances the free-energy cost of the interface between coexisting phases with the gain 
from the bulk ordering. Mathematically, the barrier is calculated by first solving a saddle-point 
or mean-field equation with the boundary conditions that the system is disordered at infinity and 
ordered at the origin, say. Now for fluctuation driven first- order transitions quite generally [13], 
there is no phase coexistence in mean-field theory so this first step of the calculation cannot be 
taken. 

Our work addresses this issue by developing a coarse-graining procedure whereby the modes 
with wavevectors in the range 

< q < qo — A, and qo + A < q < oo (1.1a) 

are averaged over, and the remaining modes with 

\q-q \<A (1.1b) 

are retained as fluctuating degrees of freedom. When A — > all modes are averaged over and we 
recover the bulk thermodynamic averages calculated by Brazovskii [1]. For A > 0, we obtain an 
effective free energy Tk for modes with wavevectors in the band (1.1b), whose average wavevector is 
qo but whose envelope can vary on a scale L > A -1 . For A not too large the bulk phases described 
by J- \ do show phase coexistence in mean-field theory, so a critical barrier height can be estimated 
from classical nucleation theory for given A. We are then left with the problem of determining the 
proper value of the cutoff A, and for this we propose a self-consistent procedure whereby A is equal 
to the local rate of variation of the envelope in the critical droplet solution. 

In cases when the ordered state is spatially anisotropic it is important to optimize the shape of 
the critical droplet, which will not be spherical due to the anisotropy of the interface free energy. 
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The standard way to carry out this optimization is known as a Wulff construction [14] and it 
arises in the present problem as well. The critical droplet we are led to for a lamellar ordered 
state has the form of a long needle with lamellae transverse to the needle axis. In addition, just 
as for ordering in liquid crystals [15], we find that under some circumstances line defects can be 
introduced into the structure and lead to deformations of the lamellae which lower the critical 
barrier height. Since these effects depend on gaining surface free energy at the cost of line energy, 
defects are only favorable very near coexistence where the critical droplet has large dimensions and 
the barrier height itself is large. 

In the Brazovskii theory there is a unique small parameter A representing the smallness of the 
dimensionless noise strength, or equivalently of the dimensionless coupling constant. The reduced 
temperature f = f c < at which the first-order transition takes place is of order |f c | ~ 

A 2/3 

below the mean-field transition which takes place at f = 0. For a d-dimensional system the critical 
free-energy barrier for creation of an anisotropic (Wulff) droplet is found to be 



B 



A (l-d)/6 



Tr 



-, d-l 



\ T \ ~ \ T c\ 
B W ~ A (5-<*)/2| f |(d~7)/ 2) 



In contrast an isotropic droplet has a barrier 



B, 



B, 



A (i-d)/3 



T > Tr 



-1 



\ T c\, 

A 2 / 3 . 



A (3-d)|-|d-4 



T \ ~ \ T c\ 

If | > If I 



(1.2a) 
(1.2b) 

(1.3a) 
(1.3b) 



from which we see that the anisotropic barrier is always less than the isotropic one. An important 
conclusion of (1.2a) is that for |f | ~ |f c | ~ A 2 / 3 the dimensionless barrier [Byy ~ A^ 1- ^/ 6 ] is still 
large for A C 1 and d > 2, i.e. the probability of nucleation is low. On the other hand this 
is the point at which the droplet size becomes of the same order as the interface width, and for 
lower quenches (|f | > |f c |) the droplet will be ramified [16]. We thus do not expect spontaneous 
nucleation of well-defined Wulff droplets. 

Although our model and methods follow those of Fredrickson and Binder [10] our conclusions 



are different. These authors only estimated the barrier for isotropic droplets in the domain ( 1.1a ), 
and did not consider either ramified or anisotropic droplets, which according to our estimates are 
the ones which are likely to be nucleated. 

In Sec. II the coarse-grained free energy is obtained, first by a phenomenological argument 
which follows the derivation of the expanded Brazovskii free energy in Ref. [10], and then by using 
a novel momentum-shell renormalization group. In Sec. Ill interface and droplet free energies are 
estimated, based first on the coarse-grained free energy with fixed cutoff, and then on a free energy 
with self-consistently determined cutoff. The self-consistency is shown not to affect the scaling of the 
nucleation barrier. The anisotropy of the interface free energy on the other hand is important, and 
it leads to an anisotropic critical droplet whose shape is determined by a Wulff construction. The 
contributions to the barrier height of defects and distortions of the order inside the critical droplet 
are estimated and it is shown that these effects are negligible except very close to coexistence, when 
the nucleation barrier and critical droplet dimensions are very large. Section III concludes with a 
brief discussion of experiment and of ways to pursue the theory in more quantitative directions. 
The derivations of the coarse-grained free energies are described in the Appendixes. 
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II The Coarse-grained Brazovskii Model 



A Bulk behavior 

Our starting point is a phenomenological model with relaxational (Model A) dynamics and a 
Brazovskii free energy which we write as 

dd = -JT + V, (2-la) 

t = fd^if^+^+^+m 2 ], (2.1b) 

(r?(x,t)r?(x',t')) = 2 5 (d) (x - 5t')S(i - F) . (2.1c) 

In the above equations the quantities qo, £o are considered to be of order unity and there is a single 
small parameter 

A <C 1 , (2.2) 
and a control parameter f. The static (long-time) solution has the bulk free energy 

$(f, A) = -ln [< e -^>] , (2.3) 

where the angular brackets in ( |2.3| ) denote an average over the Gaussian noise fj, which can be 
represented by a functional integral [see below; the effective temperature, or noise strength, has 
been scaled to unity in Eq. ( |2.1c| )1. The relation between our model and various physical systems 
can be recovered by referring to Eq. (2) of Brazovskii [1], Eqs. ( [2. 3D and (4.1) of Fredrickson and 
Binder [10], and Eqs. (A26-29) of Hohenberg and Swift [5]. The essential point is that the small 
coupling constant A reflects the smallness of the noise strength in the original systems. 

According to the derivation of Brazovskii [1], for small A the self-consistent propagator in the 
disordered phase is obtained from the Hartree diagram in Fig. la, as 

9- 1 (q)=f + $(q 2 -cf ) 2 , (2.4) 

with 

A f d d q 1 

2 J {2^Yf + e {q 2 -qlf ' 



As shown below, the solution f(f ) of Eq. Q2.5D remains positive for all f, so the linear instability 
of the disordered state, signalled by the vanishing of r in mean-field theory (A = 0), has been 
completely eliminated in the Hartree approximation. For f < there is a competing ordered 
solution with [1] 

$(x)) = ie^°' s + cc , (2.6) 



a propagator leading to 



fA=f+ ^i ( d X- -2, + a i^i 2 ' ( 2 - 7 ) 

2 J (2ir) d r A + ^(q 2 - q 2 Q ) 



and a field h conjugate to ip, given by 

1 



h = A(r A -^X\A\ 2 ). (21 
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The bulk free energy difference per unit volume A$ between the disordered and the ordered states 
can be obtained from the relation [1] 



A$ 



dA' dA 



2hdA' 



fA jdA' , 
2h—dr. 
dr' 



For small r, ta Eqs. (|2.5|) and ((2? 



where 



may be rewritten in the form 

f = t + aA/vf, 

f + a\j \J¥a~ + A|A| 2 , 



(2.9) 

(2.10) 
(2.11) 



a = (qt 2 S d ir)/4(2ir) d Z 2 = vra/2 , (2.12) 
and Sd is the surface area of the d-dimensional unit sphere. The essential point is that in all d > 2, 



the integral in Eqs. ( |2.5| ) and (2/7) contributes only near the surface of the d-dimensional sphere of 
radius %, and it is a one- dimensional integral in the radial direction. The transverse dimensions 
only contribute to the coefficient a. As shown by Brazovskii, the free energy difference per unit 
volume A«3> of Eq. ( |2.9[ ) then takes the form 



r 

2A 



or 



1/2 + ! + -f + ^L4| 4 



(2.13) 



where r(r) and r^(r) are given by (plO|) and ( gH| ) [as mentioned in HS [5] some numerical 
coefficients which were incorrect in Eq. (14) of [1] have been changed]. The free energy difference 
changes sign at 

f = f c = -2.03(aA) 2/3 , (2.14) 

which is the bulk (first-order) transition point in the Hartree approximation. Let us introduce the 
reduced variables 

(aA)- 2 / 3 f, 



r 

T 

A 
x 



(«A) 



-2/3; 



A 1 /2(«A)- 1 /3^, 
A e iq °- x + cc , 
A 1 / 2 (aA)" 1 / 3 A, 
5L/I, 

got, 

I = (2| 2 g )(aA)- 1 / 3 ^Co(aA)- 1 / 3 , 
in terms of which Eqs. (2.1b), j2l|) , (|?7|), ( fTPI and (gig) become 



/3 
J - 

r 

A$ 



^(aA) 4 /^- 1 , 

/ A {^ 2 + ^v 4 + ^r'Kv 2 + <&) 

r + ^/(2r 1 / 2 ), 
r + vr/(2r 1 4 /2 ) 



A$ 



A) 4 / 3 
A 



-2.03 (vr/2) 



--L4I 4 
A \A\ , 

(ftA) 4 / 3 



2\ , 71 r 1/2 



2/3 



-2.74. 



(2.15a) 
(2.15b) 
(2.15c) 
(2.15d) 
(2.15e) 
(2.15f) 
(2.15g) 
(2.15h) 

(2.16) 
(2.17) 

(2.18) 

(2.19) 
(2.20) 

(2.21) 
(2.22) 
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In the scaling of Eq. (2.15) the coupling constant in ( 2.18Q is of order unity and the only small 
parameter appears in the coefficient of the gradient, q^" ~ A 2 / 3 . 

It is instructive to expand the free energy difference A<£ in the order parameter \A\ 2 , retain- 
ing only the first three terms. The Hartree result (2.21) can be considered as a function of the 
independent variables r and \A\ 2 via Eqs. (2.19-20) and expanded in the form 

A$ = r\A\ 2 + ^u\A\ 4 + -^l^l 6 + 0(|^| 8 ). (2.23) 
with r = r(r) given by the solution of Eq. fl2.19Q , and 

u = (l-vr/4r 3 / 2 )/(l+^/4r 3 / 2 ), (2.24) 
w = (9vr/4r 5 / 2 )/(l+^/4r 3 / 2 ) 3 . (2.25) 

These functions are plotted in Fig. 2, from which it is seen that u becomes negative for r < 0, thus 
creating the first-order transition. 

The approximate free energy difference ( 2.23; ) vanishes for 

r = r c = 9u 2 /16w, (2.26) 

which occurs for 

t = t c = -2.51, (2.27) 

which is close to the "exact" Hartree value r c = —2.74 in Eq. ( 2.22] ). 

Equation ( |2.23| ) is precisely the same (when allowance is made for slightly different scalings) 
as Eq. ( |2,22| ) of Fredrickson and Binder [10], which was arrived at by using a "Hartree potential" 
r# [(■;/>)] in place of the "bare" free-energy functional ( 2.1 8| ) . While we believe that the proper 
physical interpretation of this potential is in terms of a partially coarse-grained potential J-^{ip) (as 
explained in the next section), for the purpose of calculating the bulk thermodynamic properties 
the coarse-graining can be carried out to arbitrarily long wavelengths (A = 0), and the potential 
A$ of Eq. (PE23D agrees with T H [(ip)] of Ref. [10]. 

Let us inquire into the domain of validity of the Hartree approximation ( [2.21 ). As mentioned 
by Brazovskii [1] and by Swift and Hohenberg [4], a simple estimate is obtained by finding the 
parameter region where the correction term in the self-energy [Fig. lb] becomes of the same order 
as the terms retained in Eq. fl2.5|) [Fig. la]. Let us write Eq. ( |2.5[ ) generally as 



f = r + S, (2.28) 

E = £h + £ 2 + 0(A 3 ), ( 2 - 29 ) 
with Tifj ~ A/f 1 / 2 the term retained in Eq. ([^), and 

S 2 ~ A 2 /f 3 / 2 , (2.30) 

the contribution from the diagram in Fig. lb. In order for the perturbation expansion in A to be 
valid we require 

S 2 ~ A 2 /|rf/ 2 < r, (2.31) 

leading to r > A 4 / 5 . In this domain the Hartree f, given by Eq. ( |2.10D , scales as f ~ (A/|f|) 2 (i.e. 
f -C |f|), and Eq. ( 2.31[) implies 

A 3/5 . (2.32) 



T\ < TGI 
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We use the subscript G in ( p.32|) in analogy to the Ginzburg criterion for validity of mean-field 
theory in critical phenomena [18]. Here, since we are dealing with a first- order transition the theory 
is self-consistent so long as 

\f c \ <C |f G i|, (2.33) 

which holds for 

A 1 / 15 < 1 . (2.34) 

We thus see that the strict self-consistency only holds for impractically small coupling strengths. 
We may also note that a less stringent criterion was invoked by Brazovskii [1], namely 



S 2 < Si = E H , (2.35) 



which replaces Eq. ( |2.32D with 



r < \ T G2\ 



A 1/2 , (2.36) 



and Eq. (|2.34Q with 



A 1/6 < 1. (2.37) 
B Nonuniform systems 

In order to study phase competition and nucleation we must be able to describe situations where 



the envelope A of the order parameter (ip) of Eq. ( 2.15d ) can vary in space, as for instance at the 
interface between ordered and disordered domains. For this purpose we introduce a coarse- graining 
procedure, whereby modes with wavevectors larger than some cutoff [19] A are averaged over using 
the Hartree approximation of Brazovskii, and modes with wavevector less than A are retained as 
fluctuating modes in averages such as Eq. (|2.3|) . In particular, let us suppose that the starting free 



energy ( 2.18j ) involves modes with wavevectors in the range < q < 2Ao = 2qo, i.e. 

< \q - go | < A , (2.38) 



with Ao = qo- The average in Eq. ( |2.3D is then a functional integral 

exp[-/3$(r)] = / V[%] exp{-/Wty,]}, (2.39) 

J[0,Ao] 

where the symbol below the integral indicates the inclusion of all modes ip(q) with wavevectors in 



the range (2.38) with upper cutoff Aq. We now define a coarse-grained free energy obtained by 



eliminating the modes in the slice A < \q — go I < Ao- This corresponds to setting 

exp[-/?$(r)] = / V[%] exp{-/3.F A [^]}, (2.40) 



with 

exp[-/3^ A ] = / V[il> q ) exp{-A7^ g ]}. ( 2 .41) 

J[A,A ] 

With these definitions we show below that 

$ = ^A=o, (2.42) 

and 

^ = ^A=A - (2.43) 
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Moreover, in Appendix B we show that if the integrals implied by Eq. ( 2.4CQ are carried out using 
the Hartree approximation of Brazovskii and if the resulting !F\ is expanded in tp 2 , then we obtain 

= J [oa] d d * [\r(A)^ + ^u(A)^ 4 + ^w(A)^ + (^) [( V 2 + qlW) , (2-44) 

where the symbol below the integral again indicates that ip{x) has variations involving wavevectors 
in the range < \q — go I < A. The coefficients are given by [19] 

r(A,r) = r + 0i(A,r), (2.45) 
u(A,r) = [l-fa(A,r)]/[l + fa(A,r)}, (2.46) 
w(A,t) = 120 3 (A,r)/[l+0 2 (A,r)] 3 , (2.47) 



with [20] 



dk 

I- [r(A.r) + (2A8) 



These quantities are evaluated by first solving the transcendental equation ( 2.45| ) for r(A, r) and 
then inserting the result into fa and fa to obtain u(A, r) and w(A,t). 

The derivation of Appendix B and the result in Eqs. (2.45-48) are closely related to those of 
Fredrickson and Binder [10], but the physical content is rather different. We are separating out 
the short-wavelength modes involved in generating the first-order transition and the possibility of 
phase coexistence [u(A) < 0], from the long- wavelength modes involved in building interfaces and 
other large distortions of the order parameter. The "Hartree potential" Tjj of Fredrickson and 
Binder is essentially J-\=o = 3>, which is no longer a functional of a fluctuating order parameter 
ip(q) where \q — go I < A. This potential depends only on the average order parameter {ip){q = qo), 
which is only a function of r in a uniform system. It thus seems to us that T\ rather than Th is the 
proper starting point for the evaluation of droplet free energies and metastability lifetimes, though 
as shown in the next section our actual results do not differ significantly from those obtained using 
Th- To make a quantitative estimate of the droplet free energy by means of ^[ip], Eq. ( 2.44| ), we 
must specify the value of A and this is done in Sec. IIIB. 

We should also mention that in principle the coarse-graining procedure applies to the full 
dynamics of Eq. ( |2.1a| ), not just to the static averages such as Eq. (|2.3|). Since our treatment 
of metastability does not go beyond the estimation of "energy" barriers [21], we have not pursued 
this question here. Clearly, however, a more complete theory would have to take into account the 
effects of coarse-graining on the dynamics and on the "entropic" corrections to the barriers height 
coming from fluctuations about the saddle point [11, 12]. 

C Properties of the coarse-grained free energy 

The model (2.44-48) defines a free energy !F\ which was designed to interpolate between the bare 
free energy T of Eq. ( 2.18; ) for [20] A = Aq = oo, and the bulk average free energy A<E> of Eq. ( 2.23 ) 



for A = 0. Indeed, it follows from Eq. ( |2.48|) that for A = Ao = oo, n (oo) = 0, so Fa=oo just 
reproduces the bare free energy T of Eq. (|2.18[) . On the other hand from Eq. (|2.48] ) f° r A = we 
have 4>i (0) = n/2r 1 / 2 (0), fa(0) = vr/4r 3 / 2 (0), fa = 3vr/16r 5 / 2 (0), so r(0), tt(0) and w(0) agree with 
the quantities r, u, w, respectively, defined in Eqs. (|2.19| ), ([2.24D and ( p.25|) . Moreover, for A = 



there is only one mode left in the free energy in Eq. (|2.44), namely the average order parameter, 



for which we make the Ansatz ( |2.15d| ). The free energy per unit volume then becomes precisely 



the expanded Hartree expression given in Eq. (2.23). 
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For intermediate A values we may evaluate r(A, r), u(A, r), w(A, r) numerically from Eqs. (2.45- 
48). The theory is of physical interest for values of A which are low enough so that r(A, r) > 0, 
implying that the disordered phase is metastable. In that case it turns out that to good accuracy 
we can represent the coefficients in the form 

r(A,r) ~ r(0,r)(l + At) = r(l + At), (2.49a) 
u(A,t) ~ u(0,t) = u, (2.49b) 
w(A,t) ~ w(0,t)=w, (2.49c) 



results which are valid for 



< A < -r~\ r < 0. (2.50) 



D A momentum-shell renormalization group 

It is natural to ask whether the coarse-grained free energy obtained in the previous section could 
not be derived in a more standard way, using Wilson's momentum-shell renormalization group [22, 
23], for example. It turns out that because the ordering in the Brazovskii model (2.1b) involves 
condensation onto the surface of a sphere in reciprocal space, the usual methods are difficult to 
implement, and various authors have found it necessary to introduce modifications of the model 
in order to obtain renormalization group recursion relations [24, 25]. However, recent work on 
the renormalization group for Fermi liquids at low temperatures, where the wavevectors of the 
important modes also lie on a sphere, suggests that a direct perturbation expansion might work. 
Making use of some of the techniques developed by Shankar [26] for the Fermi liquid we have been 
able to derive recursion relations for the Brazovskii model (2.1b), keeping essentially the same type 
of Hartree diagrams as in earlier work. The derivation is summarized in Appendix A, and the result 



is again of the form (|2.44| ) , where now r(A), u(A) and w(A) are defined by differential recursion 
relations which are quoted in Eqs. (A. 24). Solutions of these equations yield coefficients which 
are close to those of the simple approximation (2.45-48) for r > 0, and show a similar dependence 
on A and r for r < 0. There is, however, an important difference, in that the solutions of the 
differential recursion relations (A. 24) are not defined for all A and r due to singularities for r < 0, 
where r(A, r) + A 2 = 0. In particular the quantity r(0, r) vanishes at a finite t% < r c , unlike the 
solution of ( p. 45 ) which remains positive for all r. These singularities make it difficult to use the 



renormalization group to estimate droplet free energies for sufficiently negative r, so we shall rather 
use the phenomenological coefficients (2.45-48), which are defined for all A and r. Nevertheless, 
the recursion relations are well behaved for larger r (including part of the metastable range), and 
they are of intrinsic interest, so we have presented their derivation in Appendix A. 

Ill Droplet Theory of Nucleation 

A Phenomenological theory 

We shall be interested in estimating the free-energy barrier [21] for nucleation of a critical droplet 



of the ordered phase ( 2.15d ) embedded in the disordered phase ip = for r < r c . Having obtained 
an effective free-energy T\ with u < 0, we first estimate the barrier height for a critical droplet 
with fixed A, deferring to the next subsection the question of the proper choice of A. As is done in 
standard nucleation theory [11-13], we seek a saddle-point solution for J^\, i.e. a localized solution 
of the differential equation 

^y( V 2 + gg)V + r(A)V + ^u(A)V 3 + ^w(A)^ = 0, (3.1) 
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with ip = for |x| — > oo and tp ^ in the interior. In this equation the coefficients r, u, w are O(l) 
and are given by Eqs. (2.46-49), and qo ~ A -1 / 3 » 1. Given a solution ij) s of the above saddle-point 
equation, the classical nucleation barrier [21] is 

B = Ba = «]. (3.2) 
Although we have now reduced the calculation to a straightforward problem in differential 



equations, which could certainly be attacked numerically, Eq. (3.1) is still quite difficult to solve 
accurately, and we shall rather attempt to estimate the value of B\ analytically, in particular 
its dependence on the parameters, which are qo, r, u and w [in the remainder of this subsection 
we consider these parameters to be basic, and suppress their dependence on A, A and r]. The 
difficulty in finding a general solution of Eq. ( |3.1| ) arises from the existence of three widely different 
characteristic lengths. The first is the wavelength g^ 1 of the ordered phase, which in the present 
scaling (2.15) is vanishingly small. Then there is the thickness of the interface between the ordered 
and disordered phases, which is of order unity (see below, however), and finally the droplet size, 
which diverges at coexistence where the bulk free energies of the two phases are equal. 

1 Interfaces 

We begin by considering only ordered states of the form 

i>{x,y) = A(x,y)e iqoX + cc, (3.3) 

describing small distortions of a pattern consisting of parallel planes (or lines in two dimensions), 
with wavevector in the x-direction, say. Then for solutions whose envelope A(x,y) varies slowly in 
space the latter can be shown to satisfy the well-known amplitude equation [27] 

[d x + (i/2q )d 2 y ] 2 A(x, y) = 2rA + u\A\ 2 A + |^|A| 4 A, (3.4) 

where y here denotes the (d — l)-dimensional vector transverse to the direction x of alignment of 
the structure. For bulk states, with A =constant, and for u < — (4rw/3) 1//2 , the equation has two 
types of solutions, disordered [A = 0] and ordered LA = ±Aq] with 

Al = 3[(u 2 - 4ru>/3) 1/2 - u]/w. (3.5) 



Their free energies become equal when r = r c , with r c given by ( j2.26|) , though now r, n, w denote 
r(A), n(A), w(A) rather than r(0), n(0), w(0), as in that equation. 

Let us now consider interfaces between the ordered and disordered states. Clearly there are two 
types of simple interfaces, longitudinal and transverse. A longitudinal interface is a solution of Eq. 



(3.4) in which A only depends on x with A = for x — > — oo, and A = A$ for x — ► +oo say (see 
Fig. 3). This is a quintic Ginzburg-Landau equation as discussed, for example, by Fredrickson and 
Binder [10]; the interface thickness is of order 

£|| ~ (2r)- 1/2 , (3.6a) 
and the interface free-energy (per unit area) is of order 

f||~£||/o, ( 3 - 6b ) 



n 



where /o is related to the bulk free energy difference (per unit volume) 



A$ = -f S, (3.7) 
S = (r c -r)/r c , (3.8) 
fo ~ 81|m| 3 /32w 2 . (3.9) 



[If fo is defined by Eq. ( p79| ) then ( |3.7| ) is strictly correct only for 5 — > 0, but it is quite accurate 
up to 5 = 1, and we shall use Eqs. (3.7-9) for all 5 in our estimates.] 
A transverse interface, on the other hand (Fig. 4), is a solution of 

= -Ld$A(y) + 2rA + u\A\ 2 A + -w\A\ 4 A, (3.10) 
4(/q y 6 

which goes to Aq for y — * +oo and to zero for y — > — do, say. The interface thickness is now of order 

U ~ [S^r]- 1 / 4 ~ A 1 / 6 ^ 1 / 4 , (3.11) 
(which is small at fixed r for Ac 1), and the interface energy is 

o-j_~£l/o<<7||. (3-12) 

It is this anisotropy in the interface energies which leads to a preference for nonspherical droplets, 
as detailed in the next subsections (this was missed by Fredrickson and Binder [10]). 

2 Critical droplets: isotropic case 

Let us now estimate the free energy of a critical droplet solution of Eq. (|3.1|). Near coexistence (i.e. 
for 5 <C 1), the dimensions of the droplet are very large, and we can use the standard separation [12] 
of T into bulk and surface contributions to estimate its free energy. We first consider an isotropic 
droplet, i.e. one made up of concentric lamellae (we refer to three-dimensional structures in our 
discussion, but the results are applicable in d = 2 also). Then the edge of the droplet is a purely 
longitudinal interface [see Fig. 5a], and the energy of a droplet of radius R can be written in the 
standard way as 

AF{R) ~ R^a\\ - R d f S, (3.13) 
The critical droplet is the one with maximum free energy (as a function of R) , namely it has radius 

Riso~Z\\/S, (3.14) 

and corresponds to a free energy barrier 

B iso = Cf/o/^ 1 ~ [r- d / 2 \u\ 3 /w 2 }/5 d - 1 . (3.15) 

As usual the radius of the critical droplet and the barrier height diverge at coexistence, r — > r c , 
5 — > 0. On the other hand, as the undercooling 5 grows the radius decreases and when 5 = 0(1) 
the droplet size R is of the same order as the interface thickness £m, and the above estimate base 
on a separation between bulk and surface free energies becomes questionable. As pointed out by 
Unger and Klein [16] in an analogous situation, the critical barrier height may still be large at that 
point and it is useful to extend the calculation to the "ramified" case 5 = 0(1). 



Let us first recast the above calculation in terms of the free energy J^\, Eq. ( 2.44 ) and the 



associated Euler-Lagrange equation (pA\). For an isotropic configuration made up of concentric 
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lamellae the wavevector of the ordered state is radial, so the order parameter may be taken in the 
form 

ip = A{p) e iqop + cc, (3.16) 

with real A, where p is the radial coordinate (we neglect all transverse variation for this simple 
estimate). The equation for the droplet takes the form 

- 28 2 p A - -d p A + 2rA + uA 3 + ^A 5 = 0, (3.17) 

in leading order where terms neglected are of order (pqo)~ 1 d 2 A, qQ 1 d 3 A, or higher. The corre- 
sponding free energy is to leading order (in three dimensions) 

F iso iA]=S 3 J ' dp p 2 l(d p A + p^A) 2 + rA 2 + jA 4 + ^A e \ . (3.18) 
Let us make a change of variables 

A = {r/\u\) l l 2 A, (3.19a) 

p = r~ l/2 p, (3.19b) 

whereby ( 3.17| ) and ( |3.18| ) become, for u < 

2djA + {4/ p)dpA - 2A + A 3 -(3/32){l- 5) A 5 = 0, (3.20) 

Fiso = S d r^l 2 \u\- l ^ iso [A\, (3.21) 



and for d = 3 



F iso [A] = Jdpp 2 {(d~ p A + p^A) 2 + A 2 - - A A A + (1/64)(1 - <5)I 6 | . (3.22) 



where Eq. ( p.8[) was used. [An analogous formula for T can be derived for arbitrary d\. 

We now seek a (critical droplet) solution of Eq. ( 3.20Q which vanishes for p — ► oo and is nonzero 
for p = 0, with dpA = for p = 0. The only parameter left in Eqs. ( |3.20| ) and ( 3.22| ) is 5, and for 
S <C 1 the solution reaches the bulk value A 2 = A\ = 8 [see Eq. (|3.5| )1 at p = 0. The position R of 
the interface results from a balance between the surface and bulk terms in Eq. (|3.22j ), and an 
argument analogous to the one in Eq. ( 3.13 ) leads to the results TZ ~ 5~ 1 , Bi so = Ti so [A] ~ 5~ 2 
for 6 < 1. From Eq. (|3T2l] ) we then see that 

B iso = r^! 2 \u\- l B iso ~ r^^M- 1 ^ 1 ), (3.23) 

which agrees with Eq. ( |3l5| ) since for 8 <C 1, r ~ r c ~ |u| /«;. 

Let us now consider the case 5 = 0(1) when R = 0(1) and A(p = 0) no longer reaches the bulk 
value Aq(5) corresponding to (|3.5|). There is still expected to be a unique solution of Eq. ( |3~l8|) 
satisfying the boundary conditions, and since all coefficients in ( p.18 ) are of order unity, the barrier 
is given by Bi so = .Fj SO L4] = 0(1). This corresponds to the "ramified" droplet of Unger and Klein 
[16], for which the separation between bulk and surface is not applicable, but a well-defined barrier 
can still be calculated from Eq. ( p. 20 ). Thus for arbitrary 5 we have 

B iso ~ r^l 2 \u\- l B iso {8), (3.24) 

with Bi so {8) an 0(1) function proportional to 8 1 ^ d for 8 <C 1. 
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Anisotropic droplets: the Wulff construction 



All of the above estimates were based on the assumption of an isotropic solution ( p. 14 ) of the 



saddle-point equation ( |3.1| ), with concentric lamellae (or rolls in d = 2). Since, however, the cost 
a± of a transverse interface is less than that of a longitudinal interface ctii for small A [see Eq. 
( 3.12| )1, we may guess that an anisotropic solution of ( |3.1[ ) will lead to a lower barrier. The simplest 



such solution has parallel lamellae with wavevector in the x-direction, say, as in Eqs. (pj) and 



Let us first consider the classical nucleation regime 5 <C 1 and seek a solution of Eq. ( |3.1| ) 
which vanishes for x, |y| —> oo and reaches the bulk value ( |3.5[) in the center [see Fig. 5b]. Because 
of the anisotropy of the interface energy we must optimize not only the size but also the shape 
of the droplet, using the well-known Wulff construction [14]. In order to do this we consider a 
parallellipiped of length 7<T|| in the x-direction and 7<7j_ in the (d—l) transverse (y) directions, with 
7 an unknown constant to be determined (see Fig. 6). The volume of the figure is O = l d ^ d ^ 1(J \\ 
so 

7 = (n/a d ^ 1 a ll ) 1/d . (3.25) 
Then the total free energy contribution from all the surfaces is (see Fig. 6) 

T B = 2v\\(w±) d - 1 + 2(d-l)a ± ( 1 a ll )( 1 a ± ) d -\ 

= 2df rt d - 1)/d (^'t 1 ) 1/d , (3.26) 



while the bulk contribution is Tb = — /o^<5 [i n deriving Eq. ( |3.26 ), Eqs. (3.6) and ( 3.12|) were 



used]. We now maximize the total energy T = Ts + Tb with respect to the unknown volume f2, 
and find the "Wulff" values 

n = n w ~ Rfy, (3.27) 

with 

Rw ~ tw/5 ~ cr w /f 5 ~ (£„£{- V^A (3-28) 
and a total free energy barrier 

% ~ ^t'fo/S 4 - 1 ~ i d v k/5 d -\ (3.29) 
The longitudinal and transverse dimensions are 



(R\\)w ~ 7^cr||~C||A (3.30a) 

(.R±)w ~ 7VK<7± ~ (3.30b) 

( J R||)W(^±)^~% /2 ~A- 1 / 6 , (3.31) 

which diverges for A <C 1. Note also that the droplet becomes ramified (R± ~ and JF2ii ~ ^11) for 
5 ~ 1, just as in the isotropic case ( 3.14j) . 



i.e. there is an anisotropy 



In order to relate the above estimate to the fundamental equation (|3.1|), we note that for real 



d\A + [Aql)' 1 Vy A = 2rA + uA 3 + -wA 5 . (3.32) 



A Eq. (3.4) takes the form 

1 

6 

We seek a solution which vanishes for |y| — > 00, and reaches a nonzero value over a finite region 
surrounding the origin. Equation ( |3.32j ) is derived from a free energy 

T ani [A] = J dxdy I (d x A) 2 + ^2 ( V y ^) 2 + rA 2 + ]uA 4 + j^wA 6 1 . (3.33) 
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We now introduce the scaling of Eq. (3.19a) and 



x = r ~ 1/2 x, (3.34a) 
y = r- 1 / 4 (2g )- 1/2 y, (3.34b) 



which changes Eqs. ( 3.32| ) and ( |3.33| ) into 



d\A + «9|1 = 2A-A 3 + (3/32)(l - <5)i 5 , (3.35) 
Tarn = r( 7 - d ^\u\- 1 (2qo) {1 - d)/2 ^ani, (3.36) 
T ani = J dxdy |(%i) 2 + {d 2 y Af + A 2 - ii 4 + (1/64) (1 - <5)i 2 } . (3.37) 

In the classical nucleation region 5 < 1 the solution of Eq. ( |3.35| ) depends sensitively on 5, and 



the Wulff construction leading to Eq. (3.29) tells us that in the scaling of Eq. (3.34) we have 
fani ~ S 1 ^ d , Rh ~ R± ~ The solution becomes ramified (J2ii ~ R± ~ 1) when 5 = O(l) and 
in that case once again T a ni ~ 0(1). The free energy barrier ( ]3.29| ) may now be written as 

Bw ~ eil^Vo/^ 1 ~ qt d) ' 2 r- {1+d)/i \n\ i /w 2 5 d - 1 ~ ^-^V 7 ^/ 4 /!^- 1 , <5 « 1, (3.38) 

where we have again used the fact that for 6 <C 1, r ~ r c ~ |n| 2 /w. More generally we write, as in 
Eq. O, 

B W ~ qt^rV-^lur'Bwit), (3.39) 
with = 0(1), B w ~ J 1 ^ for <5 < 1. 

B Self-consistent theory 

Having estimated the barrier height for the simplest isotropic and anisotropic critical droplets using 
the phenomenological model ( p. 44 ) with given qo, r, u and w, we now address the issue of the proper 



choice of the coarse-graining scale A. The reason for introducing a coarse-grained free energy in 
the first place was that we needed to allow for variations of the order parameter on length scales of 
the order of the interface width £ and larger. Thus, in order to estimate the longitudinal surface 
free-energy an we need to impose the constraint 

A." 1 < (3-40) 

It is apparent from Fig. 7 that if q is confined to a shell of thickness 2A around the point qo%., 
then for qo 1 the transverse momentum |q y | = q y is limited by Aj_ = (2goA) 1 ^ 2 . Thus, since 
£^ = (2<7o£||) 1//2 the constraint ( |3.40| ) automatically ensures that 

AJ 1 < U, (3-41) 

which is also a necessary condition for evaluating a± consistently from In order to estimate 
bulk contributions to the droplet free energy, on the other hand, we need to let A — > as in Eq. 
( 2.23j ). These different constraints can be implemented concurrently by using a type of local-density 



or Thomas-Fermi approximation, whereby A is adjusted self- consistently to be equal to the local 
longitudinal rate of variation of the envelope A(x,y): 

A = A A (x) = A~ l d x A. (3.42) 
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We thus seek an extremum of the free energy functional 



KU^H = J dxdy {(d x Af + [^j ( V 2 y A) 2 + +r{K)A 2 + \u(A)A* + ^(A)A 6 } , (3.43) 

with A determined self-consistently at each point via Eq. (|3.42 ). It can be shown that the Euler- 
Lagrange equation is no longer precisely given by Eq. (]3.4|), but rather by 



d 2 x A - {AqlY 1 S7yA = 2f(A)A + u(A)A 3 + ^w(A)A 5 , (3.44) 

6 

with 

r(A) = r(A) + Adr/dA, (3.45a) 

u(A) = u(A) + Adu/dA, (3.45b) 

t&(A) = u>(A) + Adw/dA. (3.45c) 

Referring to Eqs. (2.49) we see that du/dA and dw/dA are small, while dr/dA ~ r so we set 

f(A) = r(0)[l + 2Ar], (3.46a) 
u(A) = u(A) = u, (3.46b) 
w(A) = w(A) = w. (3.46c) 

The self-consistent droplet is given by the solution of Eq. ( |3.44 ) , subject to the constraint ( |3.42| ) 
or equivalently by the saddle point of ( 3.43j ), subject to the same constraint. For the isotropic 
droplet Eq. Q3.42 ) is replaced by 

A A (p) = A-%A, (3.47) 

and the gradient terms i n (|3.43| ) an d (|3.44| ) are replaced by [d p A + (2/p)A] 2 and \d 2 p A + (2/p)d p A], 
respectively, as in Eqs. (13.171) and (|3.1g|) . 

We have not solved these equations numerically, but we may once again use dimensional argu- 
ments to estimate the effect of the self-consistency on the barrier height. In the isotropic case and 
for 5 <C 1 we argue that the main change compared to Eq. ( p,15[) is to select A = £7 in evaluating 
the surface term, whereas we may take A = R^ 1 ~ 5 ~ in the bulk term. Thus the barrier is 
changed to 

B% ~ti\\(Hq 1 ))fom)/S d - 1 - (3-48) 
Now from Eq. ( |3.46a ) and the relation £||(r) ~ f -1 / 2 we may show that 

£||(r(£|f))~3.2£||(f(0)) = 3.2C|| , (3.49) 

i.e. the self-consistency increases the barrier by a factor of order unity. More generally, for arbitrary 
5, we introduce the scaling (3.19), but where now r and u stand for r(A = 0), u(A = 0). Then Eqs. 
720| ) and ( |3.22j ) become, respectively 

2djA + (2/p)d rho A - 2f(A)i + A 3 - (3/32)(l - 5)A 5 = 0, (3.50) 
K C [A] = Jdpp 2 {(0 p -i + p- l A) 2 +~f (A)A 2 - \tf + (1/64)(1 - b)A 2 \ j , (3.51) 



and 
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with 



5 = [ r (A = 0)-r c ]/r c = (r-r c )/r c , (3.52) 
f(A) = l + Ar^V, (3.53) 
f(A) = l + 2Ar 1 /V. (3.54) 

The self-consistency relation ( |3.47| ) becomes 

A = A A (p) = A~ l dpA. (3.55) 

Once again all coefficients in ( 3.5C| ) are of order unity for 5 = 0(1), so Eq. ( [3,23 ) becomes 

B^-r^y^u]- 1 ^), (3.56) 

where r, u and 8 refer to A = 0, and Bf£ a = 0(1) is obtained by solving ( 3.5CQ subject to the 
constraint ( 3.55 ). 

In the anisotropic case the arguments go through in the same way, and they lead to the replace- 
ment of ( gjp for 5 < 1, by 

B% ~ e N [r(A = q 1 )}^- x [r(A = ^l/oKA = i^ 1 )]* 1 ^. (3.57) 
Now for 5 <C 1, -R^ 1 ~ 5 — » 0, and once again Eqs. (2.49) and ( |3.6a| ) lead to 

r(A = ^ 1 ) = 0.1r(A = 0). (3.58a) 

and 

e||(A = ^ 1 ) = 3.2^||(A = 0). (3.58b) 

Thus the self-consistent anisotropic barrier is increased with respect to ( 3.2S| ) by a factor of order 
unity, and Eq. ( |3,38j ) still holds, with r, u and 5 now referred to A = 0. 

For 5 of order unity the scaling of Eqs. ( |3.19a ) and (3.34) goes through with r = r(A = 0), and 
Eqs. ( p. 35 ) and ( 3.37 ) are replaced by 

dp + d\A = 2f(A)i-i 3 + (3/32)(l -5)A 5 , (3.59) 

J dxdy |(a £ i) 2 + (d 2 y A) 2 + ~r(K)A 2 - ^A A + (1/64) (1 - <5)i 6 } , (3.60) 



^ani [A] 

with the self-consistency relation 



A = Ar(i) = A^diA, 



(3.61) 

and r(A), f(A) given by Eqs. fl3.53|) and (|3.54|) , respectively. Then the self-consistent estimate of 
the anisotropic barrier has the same form as Eq. (3.39), 



B 



11' 



u\- L B%{5), 



(3.62) 



but now r, 5 and u refer to A = 0, and B^{5) = 0(1) is obtained by solving the self-consistent 
equations (3.59-61). 

The final results may then be expressed in terms of the original variables of Eq. (2.1b), namely 
f and A, by using the formulas of Sec. IIA. We defer a detailed examination of the results until Sec 
HID below, but it should already be clear that the anisotropic barrier ( |3.62 ) will always be smaller 
than the isotropic one ( 3.56|) , due to the factor ~ ^C^ 1 )/ 6 <g \ an d the higher power of 

r < 1 in Eq. 
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C Defects and distortions 



The calculation of the anisotropic barrier was based on the simplest ansatz designed to take ad- 
vantage of the favorable transverse interface energy a±, namely parallel lamellae with wavevector 
q ~ (/ox, as in Eq. (|3.3j). It should be clear, however, that under certain circumstances a lower- 
energy configuration can be achieved by distorting the lamellae in regions where they form (un- 
favorable) longitudinal interfaces, in order to gain additional surface free energy. In general such 
distortions tend to produce defects, whose free energy cost is linear in the size of the structure, so 
distortions are typically favored for large droplets. 

Our discussion of defects and distortions draws heavily on the work of Cross [28] and of Fournier 
and Durand [15]. In particular, the latter authors showed that for smectic liquid crystals it is 
advantageous to introduce defects in the form of focal conies, and to pack these into overall spherical 
shapes in the asymptotic limit of an infinite droplet (Fig. 8). As shown below, the same arguments 
hold for the present system but it turns out that focal conic defects are only favorable in a small 
region near coexistence (5 <C 1, t — * t c , R — * oo). Similarly, we have considered overall distortions 
of the anisotropic droplet into an annulus, in order to eliminate the costly longitudinal interfaces 
at the tips (see Fig. 9). This introduces dislocations in the bulk of the structure, however, and 
according to our estimates it is not favorable for any value of r. 

Unless otherwise noted we will consider three-dimensional systems in this section, though similar 
arguments can be given for d = 2, or any other d > 3. Following Cross [28] let us write the gradient 
term in the free energy fl2.18|) in the form 



8?o 



where 



J d 3 x[(v 2 + ^] 2 = ^/rf 3 XK[(vft) 2 +4(^) 2 ], (3.63) 



^ = Ae i( + cc, (3.64) 
VC = (qo + Sq)h, (3.65) 



and the bending constant is given by 



K = ^A 2 (3.66) 
in these units. 

Let us now consider the effect of inserting a focal conic into the Wulff droplet, as sketched in 
Fig. (8a,b), and let us denote by L the average radius of curvature of the bend in the structure. 
This length also corresponds to the length of the core of the disclination defect line, which has an 
average cross section £ 2 . Then as argued by Fournier and Durand [15], the introduction of the focal 
conic yields a bulk energy cost of bending of order nL~ 2 per unit volume (or kL overall), and a 
term «Zln(L/£|i) to account for the core of the defect. The surface energy, on the other hand, is 
decreased by an amount ~ AaL 2 where 

Act = a FC — aw (3.67) 

is the difference between the structures in Figs. 8b and 5b, respectively. Now according to the 
calculation in Sec. IIIA3 we have aw ~ /o^y^i.) 1 ^ 3 ) whereas for the structure with a focal conic 
most of the interface is transverse, so dpc ~ fo€± *C The free energy difference between the 
two structures scales as 

Fw ~ Ffc ~ kL[ Ci + c 2 ln(L/M] + AaL 2 , (3.68) 
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where c\ and C2 are constants of order unity. The above quantity vanishes for a size L of order 



Lq = k/\Act\ ~ k/uw, (3.69) 

which means that it is favorable to introduce focal conies into the Wulff droplets for Rw = 
&w/fo5 > L , i.e. for 

X ^ X ° w a W -2/3 in vn\ 

5 < 5 FC ~ T7 -f q , (3.70) 



where we have used the relations aw ~ £,w ~ (£||)£± ■> Eqs. fl3,66|) and ( 3.11 ) and the fact that in 



the present scaling and near coexistence (5 <C 1) all the coefficients r ~ r c , u, w, /o, k ~ A 2 ~ r/u 
are of order unity and only factors of 5 and qo ~ A -1 / 3 ^ 1 need to be considered. In the range 
( |3.70| ) the barrier is then determined by a±, i.e. it scales as 

Bfc ~ ~ ?o 3/2 /^ 2 , ^ < ?o 2/3 , ( 3 - 71 ) 



rather than (|3.38[) 



The above calculation assumed that the Wulff droplet was unchanged in shape and size, and 
only the texture of the lamellae was modified. Following Ref. [15], we may consider a composite 
droplet consisting of a spherical array of conical domains, each one of which contains a focal conic 
defect (see Fig. 8d). Such a shape was found to predominate for large droplets in the case of 
smectic liquid crystals [15]. The surface free energy is now still transverse and it is achieved over 
a spherical surface (a±R 2 ), since the cost is entirely in line energy (~ RlnR) and is negligible for 
large R. Thus the free energy is given by 

F*ph ~ <?±R 2 - SfoR 3 , (3.72) 

leading to 

Rsph ~ t±/8. (3.73) 
This is more favorable than the single conic of Fig. (8a, b) for R sp h > Lq, i.e. 

5 < 5 sph = a±/f L ~ a±a w ~ qo 5/6 . (3.74) 

The barrier for this droplet has the same scaling as the one with a single focal conic, but presumably 
with a smaller coefficient in the range ( |3.74| ), 

B sph ~ Qo 3/2 S~ 2 , B sph <B FC , S<q 5/6 . (3.75) 

Finally, let us estimate the free energy barrier for creation of an annular droplet (Fig. 9), which 
lowers the surface energy at the cost of splay and/or defect energy in the bulk. We first consider an 
undefected structure in which the bending of the lamellae leads to a change in the local wavevector, 
and a bulk energy cost given by the last term in Eq. ( |3.63 ). If we denote the radius by R\ and the 



width by R2 Ri, then the change in wavevector from the inner to the outer rim is 5q ~ qoR2/R\, 
and the energy cost is (5q) 2 R\R\ ~ q^R^/Ri. The surface energy is entirely transverse, so the total 
energy of the annulus is 

Fannu ~ ^O^/^l + % l ^ RlR-2 ~ RlR^, (3.76) 

— 1/2 

where near coexistence (r ~ r c , 5 <C 1) we may again set /o ~ 1, a± ~ £± ~ q . In order to find 
the critical size and shape we need to maximize ( |3.76D at fixed ratio 

ri = R 2 /Ri, (3.77) 
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and then minimize the result with respect to the ratio r/. We find r] c ~ f^ 1 / 2 ^ 1 , R\ c ~ % $ 3 ^ 2 
and a barrier 

Bannu ~ ^V 5 / 2 , (3.78) 

which is always larger than the Wulff value B\y ~ crvy/<5 2 ~ £j_/5 2 ~ l/go^ 2 , for <5 < 1. 

Thus the cost of changing the local wavevector is too high so we shall keep the constraint q ~ qo 
on average, and attempt to achieve bend by introducing dislocations in the bulk of the structure 
(Fig. 9b). Since dislocations comes form lines in three dimensions or points in two dimensions, the 
2d estimate of Cross [28] for the number Np of dislocations applies in 3d also, namely 

N D ~q R 2 . (3.79) 

Each dislocation costs an energy /o#£||£_l per unit length of core and has length R 2 , so the bulk 
energy cost of the dislocations is ~ foS£\\£±R%<Io- Thus the structure in Fig. 9b has energy 

J~d ~ Sql /2 Rl + q 1/2 RiR 2 - RiR 2 2 5. (3.80) 
The same procedure as above then yields r\ c ~ (qoS)' 1 , R\ c ~ q^ 2 \ and a barrier 

B D ~( q y 2 8)- 1 . (3.81) 

The above estimates are only valid if R 2 < R\, i.e. r/ c < 1 which means 

5 > q^ 1 . (3.82) 

In the next section we put together all the estimates to find the most favorable structure for each 
region of 5 (or r). 



D Results 



In Eqs. ( gjj ), ( ^62|) , (|3?75|) , (|3^ ) and (|^8|) we have presented estimates for the barrier heights of 
the droplets shown in Figs. (5a,b), (8b,d) and (9b), expressed in the scaled units of Eq. (2.15), as a 
function of qo ~ A -1 / 3 and the bulk parameters r, u, w, 5 = (r c — r)/r c (evaluated at A = 0), which 
depend on r. Near the bulk transition (5 <C 1) only the parameters i5 < 1 and qo ~ A -1 / 3 3> 1 
survive since r ~ r c , u and w are 0(1), and 



5 = (r c - r)/r c (\t\ 
Far below the transition (r — > 0, 5 — > 

r(r) ~ |t| -2 



-N)/|r c | = (|f|-- 
1, |t| > 1, |f| > |f c 
-(|r c |/|r|) 2 ~A 4 / 3 |7 



l)/|r c | « 1, 
we have 

" 2 , |r| > |f c 



(3.83a) 



(3.83b) 



Let us first compare the Wulff droplet with the various defected structures discussed in the 
previous section. According to Eqs. ( |3.70| ) and ( 3.74 ) the focal conic and spherical structures are 
favored over the Wulff droplet for 5 < q 2 ^ 3 ~ A 2 / 9 , and 5 < q 5 ^ 6 ~ A 5//18 , respectively. Comparing 
the barriers Bpc ~ B sp h in Eqs. ( [3.71 ) and ( 3.75 ) with the barrier for the defected annulus Bp 
in ( 3.81| ), we see that the annulus would be favorable for 5 < q^ 1 ~ A 1 / 3 , but this is precisely the 
region given by Eq. ( 3.82 ) where the annulus is no longer well defined. We conclude that it is 
only in a vanishingly small region near coexistence (when the droplet size diverges and the barrier 
heights are very large) that defects come into play, and when they do it is in the form of focal conic 
structures (see Fig. 10). 
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We are thus left with the Wulff droplet as the preferred one over most of the metastable range 
of t. In order to estimate the probability of nucleation we need the barrier By/ in the original units 
of Eq. (2.1b), since the lowest order estimate for this probability is the saddle-point contribution 
exp[— By/] (the noise in Eq. ( 2.1c|) is normalized to unity). According to Eqs. ( 2.15h ), ( 2. 16| ) , 
(|2~l7j) and (KM we have 



B 



W 



(3B W ~ \^- d y\^ d)/ \«- d ^\u\- 1 B w {5). 



(3.84) 



Noting that B w {8) ~ d 1 ^ for 5 <C 1 (r r c ), and B w {5) = O(l) for 5 = 0(1), and taking into 
account the asymptotic estimates in Eq. (3.83) we find 



By/ ~ 

By/ ~ 
By/ ~ 



d-X 



A (l-d)/6 

L|T|-|T ( 

A (i-<D/ 6j | f |/| fc | =0(1)) 

A (5-d)/2|-|( cf -7)/2 \ f \^\ f \^ 



A 2 / 3 . 



A similar calculation for the isotropic droplet of Eq. ( |3.56| ) yields 

d-l 



Bis 

Bi s 
Bis 



A (i-d)/3 



A (l-d)/3 ) 
T 



A (3-d)|^i</-i 



r l - l r c| 

\f\/\f c \ = 0(l) 

f \ > |f c |. 



(3.85a) 

(3.85b) 
(3.85c) 



(3.86a) 

(3.86b) 
(3.86c) 



From the above estimates it is clear that the anisotropic Wulff droplet has a lower barrier than 
the isotropic droplet, but the barrier nevertheless is still large [A^ 1- ^/ 6 S> 1] up to and beyond 
the region |f| ~ |f c | when the droplet becomes ramified according to Eq. (3.30). Indeed, we 
may estimate the barrier height at \f \ ~ \tgi\ ~ A 3 / 5 , where according to Eq. ( |2.32| ) the Hartree 
approximation breaks down, and we find 



Bw(f G i)~X {2 - d)/5 . 



(3.87) 



Thus for d = 3 the Wulff barrier is still large (A -1 / 5 3> 1) at this point, which in reduced units 
corresponds to the region (|r| ~ l^l/l^cl ~ A -1 / 15 3> 1). If we continue to use the asymptotic 
estimate Eq. fl3.85c|) beyond its formal range of validity we find that (for d = 3) the barrier height 
is of order one at 

|f | = |w| ~ X^W-V ~ A 1 / 2 , (3.88) 



which represents a crude estimate of the condensation point. We also note from Eq. ( 3.87] ) that 
for d = 2 the barrier height becomes of order unity precisely at f = tq\. 

Let us compare our results with those of Fredrickson and Binder [10]. As mentioned in the 
Introduction these authors confined themselves to isotropic droplets, so we should compare with 
our Eq. (3.86). The scaled free energy ( |2.10| ) of Ref [10] is the same as our T , Eq. (2.1b), since all 
quantities are O(l) except for the coupling constant 



u 



N- 1/2 ~ A < 1 



(3.89) 



[here and below A refers to our coupling constant, not to their parameter A, which is O(l); note 
also that the quantity 5fb defined in their Eq. (3.12) corresponds to our |f| — |r c | = £|f c |]. The 



estimate of the barrier in their Eq. ( jug) is B iso ~ N^ 3 [S\f c \]- 2 ~ \- 2 ^5- 2 , which agrees 
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with Eq. (3.86a) for d = 3. However, the condition of validity of the expansion in their Eq. 
fl3.14| ) is 5 F b < N~ 1/3 , not <5p B < 1 [i.e. 5 < 1, not 5 < A" 2 / 3 ]. Thus the estimate of a 
kinetic limit of metastability or condensation point where Bi so ~ 1 given by these authors, namely 
5fb ~ iV -1 / 6 ~ A 1 / 3 (or |f con d| ~ A 1 / 3 ) differs significantly from our estimate for the isotropic 
droplet |r con( i| ~ 1 which follows from Eq. ( 3,86d ). In any case, according to our calculations the 
critical droplet is an anisotropic ramified structure, and it leads to the estimate Eq. (3.88) for the 
kinetic limit of metastability. 



IV Conclusion 

In this paper we have presented estimates of the critical droplet free energy of the Brazovskii model 
(2.1) in the metastable phase of its fluctuation induced first-order transition in the weak-coupling, 
low- noise limit A < 1. Our work builds on that of Fredrickson and Binder [10], but finds a more 
favorable configuration for the critical nucleus than their isotropic droplet, by taking into account 
the anisotropy of the lamellar ordered state. Our derivation is also physically more plausible than 
that of Ref. [10] in our opinion, since the free energy we use distinguishes between the short-scale 
fluctuations leading to the first-order transition and the long-wavelength fluctuations necessary to 
build interfaces and droplets. In the final results, however, this added self-consistency changes only 
the numerical factors and not the basic scaling of the barrier height in the small parameter A. Our 
main conclusion, which differs from that of Ref. [10] is that on the scale of the first order transition 
[|^| — \t c \ ~ |f c |] the critical barrier is large, Bw ~ A^ 1 ^/ 6 » 1, for d>2. 

At this point none of our calculations can be claimed to yield realistic estimates of the lifetimes 
of metastable states in the systems described by the Brazovskii free-energy (2.1b). This is first of 
all because our results are based exclusively on dimensional analysis with no account of numerical 
coefficients, and with highly idealized asymptotic conditions, the most extreme of which is the 
one in Eq. ( |2.34| ). In addition, we have only addressed the question of finding the energy of the 
saddle-point configuration, which is the first step to estimating the lifetime of a metastable state, 
but by no means the whole answer. 

Our work is to our knowledge the first attempt to modify classical nucleation theory [11,12] 
in a controlled way to deal with metastability in fluctuation-driven first-order transitions [13]. A 
straightforward extension would be to solve the self-consistent equations fl3,59| ) ( [3.61 ) for the Wulff 



droplet numerically and thereby calculate the function B^y(5) explicitly. A more difficult (and 
more interesting) step would be to justify the heuristic arguments leading to Eqs. ( [3.42|) and (3.44) 



by a formal calculation analogous to Langer's [12] field-theoretic derivation of the lifetime of a 
metastable state in transitions where phase coexistence already appears at the mean-field level. 

The most promising application to an experimental system is to the microphase separation tran- 
sition in symmetric diblock copolymers [6-9], which was the main focus of the work of Fredrickson 
and Binder [10]. With all the caveats expressed above, we can say that we expect the disordered 
phase to be metastable against homogeneous nucleation of droplets down to a reduced temperature 
of the same order as or larger than the shift \f c \ ~ A 2 / 3 between the mean- field and the actual 
(thermodynamic) transition (|r con d| > |f c |). If the ordered phase is homogeneously nucleated we 
expect the critical droplets to be ramified [16], since the more regular Wulff needles shown in Fig. 
5b have large barriers. The actual shapes that would be observed experimentally can only be de- 
termined if one studies the subsequent evolution of the droplets as they aggregate and coarsen [10], 
a question whose elucidation requires further work. 

Another possible experimental application is to Rayleigh-Benard convection [4,5,27,28], perhaps 
near the fluid critical point where fluctuation effects are expected to be large [29]. In that case, 
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however, apart from non-Boussinesq effects [27] which might mask the fluctuation contributions, 
it may be important to take into account the effects of multiplicative noise [30] on the transition, 
since the latter could be much larger than the additive (thermal) noise. It is not clear at present 
how much of the present theory would be relevant in the convection system. 
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93ER14312. 
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Appendix A: Renormalization Group Recursion Relations 



We seek to obtain a coarse-grained free energy by constructing a momentum-shell renormalization 
group in the manner of Wilson and Kogut [22], except that we eliminate shells surrounding the 
sphere |q| = qo, rather than shells surrounding a single point in g-space, as was done in Ref. [22]. 
To do this we borrow from the techniques developed by Shankar [26] for the Fermi liquid. Rather 
than starting from the free-energy (2.1b), we generalize in the usual way [22] to include higher 
interactions, i.e. we start from 

T[i>\ =F 2 +Fa + F% + (A.l) 



* = 4! 



etc, where 



T 2 = -f dld2u 2 {l,2)(2Ti) d 5(l + 2)^(1)^(2), (A.2a) 

2 J[0,K ] 

\ [ dl...di n 4 (l, 2, 3, 4) (2ir) d 5(l + 2 + 3 + 4)^(1)^(2)^(3)^(4), (A.2b) 
i! J[0,Ao] 

I _ = tA^ / _ dMt 1 I dn u (A.3) 

J[0,Ao] ( 27r ) ^0<|?i-go|<Ao J 

and dQi is the element of solid angle in d-dimensions. We now integrate over the modes 

A /b < q- q < A , (A.4) 

and the corresponding inner shell with q < qo, and then rescale all momenta so that Ao/6 — > Ao- If 
we define 

k = q-qo, (A.5a) 

k' = bk, (A.5b) 

= b~ 3/2 ^{k), (A.5c) 

then in lowest order in ~ A we may write 

u 2 (l,2) = [f + iM-Qo)% (A.6) 

and the change in f is given by 

f' = fe 2 [f + A 2 ], (A. 7) 

where A 2 is the contribution from the Hartree diagram in Fig. la (see below). Let us now examine 
the change in U4. By an argument similar to the one given by Shankar [26] (Sec. V) we find that 
at "tree-level" the change in U4 is given by 

u' 4 (l', 2', 3', 4') = e -(6-l)9o[|fli+fl2+A3|-l]/Ao 6 3-^ 1) 2) 3; 4); ( A _g) 

where hi is the unit vector qi/|q«|. Thus by iteration of the renormalization group, only couplings 
with vectors satisfying 

|ni+n 2 + n 3 | = l (A.9) 
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will remain finite. It can be shown that for both d = 2 and d = 3 this implies that the four 
wavevectors in u±(l, 2, 3, 4) must be equal and opposite in pairs when their magnitude goes to go- 
Thus only 

u 4 (ni, -ni, n 2 , -n 2 ; q ) = u(n x ,n 2 ) (A.10) 

survives, and it satisfies 

u'(ni,n 2 ) = 6 3 [tZ(ni, n 2 ) + A 4 (n 2 , n 2 )], (A.ll) 

where A4 is the contribution from the diagrams in Figs, lc and Id. It turns out that u(fii,n 2 ) 
develops a dependence on the angle between fii and n 2 under the action of the renormalization 
group. For a general nonzero angle u is practically constant and we define 

u(ni,n 2 ) = u a . ( A -12) 
For the special case of parallel fii and n 2 (to within an angle of 0(\/F), we define 

u(hi,nx) = u b . (A.13) 

Similarly, we have 

u 6 (l, . . . ,6) -> T2 6 (ni, -ni, n 2 , -n 2 , n 3 , -n 3 ; q ) = w(&i, n 2; n 3 ), (A.14) 

^'(ni,n 2 ,n 3 ) = 6 4 [w)(ni, n 2 , n 3 ) + A 6 (ni, n 2 , n 3 )], (A.15) 

with Ag given in Figs, le and If. Under the action of the renormalization group w also develops 
an angular dependence; for general angles between ni,n 2 and n 3 we define 

w = w a . (A. 16) 
When only two of the unit vectors are parallel we define 

w = Wb, (A. 17) 
and for all three unit vectors parallel we define 

w = w c . (A.18) 

Let us now evaluate the contributions A 2 , A4, Ag, using the lowest relevant order in M4 ~ u ~ A 
for each A,-. We find 



d r o+A ° j==d-l /jn «4(n, -n,ni,-n 



A 2 = (27r)- d / dq qf- 1 / dQ 



?o+A /6 J r+^{q 2 -ql) 2 

^ a £/^' (A - 19) 



where a is given in Eq. ( ^.12 ) and £0 is defined in Eq. ( [2.15h ). In particular there is no momentum 



dependence to this order, so the second term in square brackets in Eq. ( |A.6| ) is unmodified. In 
order to find differential recursion relations we set 

b=l + t, (A.20) 
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with I <C 1, and find from Eqs. (|A.7|) and (A.l 



— = 2rH jtto- (A. 21a) 

The contributions to A4 are given in Figs, lc and Id. When we examine the contribution to u a from 
Fig. lb we find there is one channel in which both intermediate Green's functions have arguments 
in the shell being integrated over, whereas for u\, there are two. This allows the evaluation of A4 
and leads to 

dua _ „_ a£ou 2 a A | a£ w a A 

d£ ~ iUa (f + W + f + ^Ag' ( } 

and 

dub _ „_ 2g^og|Ao agoWbAo , 
b W ^lA|- (A - 21C) 

Now consider Ag. There are 15 ways of distributing the 6 arguments of uq over the external legs 
in both Figs, le and If. Examination of the number of channels which contribute to w a , Wb and w c 
leads to 

dw a 2a£ «aA 3a>£ u a w a A 

¥-** + f+W" (f + gA8)" (A - 21d) 

and two similar equations for and u; c . 

These equations have initial values corresponding to Eq. (2.1b), 

f(£ = 0) = r, (A.22a) 
u a {l = 0) = = 0) = A, (A.22b) 

w a (£ = 0) = w 6 (^ = 0) = «j c (^ = 0) = 0. (A.22c) 

Because of the scalings we have chosen in Eq. (A. 5) the quantities f(£), u ai {,(£), w a ^ c {£) all grow 
with I. To find finite quantities we make a change of variables 

A = A ^ (a\)- 1/3 e- e , (A.23a) 

r(A) = r(£)(a\)- 2 / 3 e- 2e , (A.23b) 

u a ,b( A ) = u a , h (e)\e- 3£ , (A.23c) 

w a , b , c {A) = » fl ,6, c W« 2 (aA)- 4 / 3 e- 4< , (A.23d) 



leading to 



dr u a 



dA r + A 2 ' 

du a u\ w a 



dA (r + A 2 ) 2 r + A' 

dub 2u 2 Wb 



dA (r + A 2 ) 2 r + A 2: 



(A.24a) 
(A.24b) 
(A.24c) 



2(3 



dw a ~^ u a 3u a w a 

~d~A = (r + A 2 ) 3 + (r + A 2 ) 2 ' (A.24d) 

dW b _ ~4U| UgW b ^ ^UgWg . 

dA ~ (r + A 2 ) 3 (r + A 2 ) 2 (r + A 2 ) 2 ' 1 j 

dA (r + A 2 ) 3 (r + A 2 ) 2 v ; 

with initial values at A = Ao^o^A) -1 / 3 = oo, 

r(oo) = r, u a ,6(oo) = 1, w a ,b iC (oo) = 0. (A.25a,b,c) 

Let us reexpress the free energy in terms of the running couplings u n (£) . Since the upper cutoff 
Ao is restored to its original value at each stage of the recursion (|A.7| ) , the effective free energy still 
has the form (A. 1-2) when expressed in terms of the u n (£) and the corresponding Let us now 

change to the scaling (A. 23), supplemented with the relation 

tf;(£) = v5(£)A 1/2 (QA)- 1/3 e 3 ^ 2 . (A.26) 

Then it is straightforward to show that the free energy takes precisely the form (2.44), with a 
variable upper cutoff A, and coefficients satisfying Eqs. (A. 24). The resulting functions r(A), u(A) = 
u (A),w(A) = w c (A) are similar, but not identical to the ones defined in Eqs. (2.45-48) and obtained 
phenomenologically in Appendix B. In particular, contrary to the latter functions, the solutions of 
(A. 24) do not exist for all r and A, due to singular denominators for r < t\ = —2.66. This is to be 
contrasted with the bulk transition predicted by (A. 24) at t£ 9 = —2.56. An analysis of Eqs. (A. 24) 
shows that the coefficients are undefined for 

A < A x (r), (A.27) 



where A\(t) is plotted in Fig. 11. In the domain QA.27 ) represented by the shaded region of 



Fig. 11, the quantity r(A,r) + A 2 appearing in the denominators in (A. 24) has zeroes, so the 
recursion relations cannot be solved with initial conditions at A = oo. Of course the singularities in 
(A. 24) occur for r(A, r) < 0, i.e. in a region where the disordered state is not even metastable, so 
they are not physically significant. However, in the renormalization group formulation the coarse- 
grained free energy at small A is found by integration starting from a bare free energy with large 
A, so the recursion relations are needed even in the unphysical parameter range with r < r c , A 
large, and r < 0. The phenomenological coarse-graining procedure of Eqs. (2.45-48) avoids the 
singularity, but even if it did not, this would not affect the utility of the method. This is because 
the phenomenological equations are given in integrated form, so the coefficients r(A), u(A), w(A) 
at a particular value of A are obtained by solving equations such as (2.45-48) involving only the 
same value of A, rather than by integrating down from large A. This means that the physically 
relevant small-A behavior can be obtained independently of any large-A singularities. 

In Fig. 12 we show the numerical solution of the recursion relations (A. 24) in a form similar 
to that in Fig. 2. The results only exist to the right of the shaded region in Fig. 11, i.e. for 
r > n when A = 0. The coefficients agree with the phenomenological ones for r > 0, but they 
differ quantitatively for r < due to the vanishing of r at t\. Nevertheless, for r > t\ the recursion 
relations provide a justification for the phenomenological theory, since they have been derived by 
more or less standard diagrammatic methods. 
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Appendix B: Phenomenological Coarse-graining Procedure 



In this Appendix we reformulate the derivation of Fredrickson and Binder [10] for the effective free 
energy, replacing the averages in their Eq. ( |A.3j ), which are over the whole range < \q — qo\ < Aq, 
by averages over the restricted range < \q — qo\ < A. The coefficients f(A), u(A), w(A) are just 
the derivatives 

r An = (B.i) 

for n = 2,4, 6, respectively. They are obtained from the diagrams in Fig. 1 except that the integrals 
are confined to the range < \q — qo\ < A, and no rescaling of the momentum is performed. As we 
noted in our discussion of the renormalization group in Appendix A, the number of channels which 
contribute to T4 and Tq depends on the angles between the wavevectors which are their arguments. 
We will eventually want the free energy to be a functional of a tjj of the form ( |3.3|) . Thus we want 
the wavevectors occurring in the arguments to be either parallel or antiparallel. In the notation of 
Appendix A we then want to identify with u and w c with w. 
For r the Hartree graph in Fig. la gives 

HA) = f + ci x£° * . (B.2) 



For u a the graph in Fig. lc yields 



f A ° - dk 

u a = x-u a xj K afo ( _ + @j£a)a . (B.3) 



For ut> this graph yields 

r R o dk 



u b = X-2u a Xj K *&> (f + ^ )2 - (B.4) 

The coefficient of the sixth order term, w c , is given by the graph in Fig. le which yields 

Mo _ dk 

" C = 12 A ^°FTW (R5) 

When Eqs. (B.2-B.5) are reexpressed in terms of the scaled units of Eq. (2.15) they then yield 
precisely Eqs. (2.45-48). 

Let us verify that the denominator in the (f> n , Eq ( 2.4S| ) remains positive 

r(A,r) + A 2 >0, (B.6) 

from which it follows that the coefficients r(A, r), u(A, r) and w(A, r), Eqs. (2.45-48) are well- 
defined for all A and r. The question only arises for r < and r < 0, so we set 

r = -A 2 (l - r]) 2 , (B.7) 



and ask whether rj can vanish. If we carry out the integral 0i in Eq. ( 2.48j ) we find (assuming 
T) > 0) 

which for r/ — > becomes 

In- = 2A(|r| - A 2 ). (B.9) 
V 



Since for fixed r the rhs of ( |B.9| ) is bounded by 4(|r|/3) 3 / 2 as a function of A, we conclude that 
has no solution for rj — > 0, and that the inequality (|B.6|) is always satisfied. 
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Figure Captions 



Low-order diagrams in the weak-coupling expansion of the Brazovskii model (2.1b). 
(a) Hartree diagram contributing to r, or A2. 



(b) The second-order self-energy diagram entering £2 in Eq. ( 2.29 ). 

(c) A diagram contributing to A4. The three ways of distributing the arguments of 1x4(1, 2, 3, 4) 
are shown. 

(d) Another diagram contributing to A4. 

(e) A diagram contributing to Ag. The fifteen ways of distributing the arguments of 
uq(1, 2, 3, 4, 5, 6) are not shown. 

(f) Another diagram contributing to Aq. 

2. The coefficients of the bulk free energy (|I23|), plotted as a function of r. (a) The coefficient 
r, obtained by solving Eq. fl2.19| ). (b) The coefficient u, from Eq. ( |2.24[ ). (c) The coefficient 
w, from Eq. ( |2.25 ). The bulk coefficients correspond to the limit A = of the coefficients of 



Eqs. (2.45-48). 

3. Schematic diagram of a longitudinal interface between the lamellar (x > 0) and disordered 
state (x < 0). (a) Sketch of the lamellae; (b) Order parameter \A\ vs. x. The width of the 
interface is £y. 

4. As in the preceding figure, but for a transverse interface, with lamellae for y > and a 
disordered state for y < 0. The transverse interface is thinner than the longitudinal one 

(ei<£u)- 

5. Schematic diagram of critical droplets, (a) Isotropic droplet with concentric lamellae, (b) 
The anisotropic Wulff droplet, consisting of lamellae perpendicular to x. 

6. The parallelipiped entering the Wulff construction. 

7. A portion of the Brazovskii sphere in reciprocal space, showing the relationship between the 
longitudinal cutoff A and the transverse cutoff (2q , oA) 1 / 2 3> A, for A <C qo- 

8. Schematic representation of focal conic defects, (a) and (b) represent a single focal conic 
introduced into the Wulff droplet; (c) is a biconical domain in which focal conies are juxta- 
posed; (d) is a focal conic spherical network of cones representing the equilibrium shape for 
arbitrarily large droplets. [Prom Ref. 15.] 

9. Possible distortions of the Wulff droplet to eliminate longitudinal interfaces, (a) An annular 
shape with no defects; (b) the same, but with the bend relieved by the introduction of 
dislocations. 

10. A summary of the metastability properties of the Brazovskii model in three dimensions, for 
A <C 1. Shown are both the scale of the control parameter f of Eq. (2.1b), and of the 
reduced quantity r ~ f A~ 2 / 3 . The fluctuation induced bulk transition occurs for r < 0, at 
t = t c = 0(1) [f = f c = 0(A 2 / 3 )]. The defected critical droplets are only favorable in an 
infinitesimal region near r c , 5 = [\r\ — \t c \]/\t c \ ~ A 2 / 9 -C 1. The anisotropic Wulff droplet 
is favorable for 1 > 5 > A 2 / 9 . For |r| 3> \t c \, [\f\ 3> \f c \] the Wulff droplet is ramified 
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since the interface width is of the same order as the droplet radius. The Ginzburg criterion 
beyond which the perturbation theory no longer holds occurs for |r| = |tgi| ~ A -1 / 15 3> |r c |, 
[|^l = I^Gil = ^ 5 ^ frcWi & t which point the critical barrier height is Bqi ~ A -1 / 5 3> 1. The 
barrier becomes of order unity at |r| = \T con d\ ~ A~ 1//6 , [|f| = \ f CO nd\ ~ A 1 / 2 ]. 

11. The function Ai(r) at which the recursion relations (A. 24) have a singularity. The coefficients 
are well-defined for all A if r > t\ and for A > Ai(r) if r < n. 

12. The bulk coefficients obtained from the recursion relations (A. 24). (a) The coefficient r(A = 
0,t) vs. t; (b) u(t) = iift(A = 0, r) vs. r; (c) io(r) = w c (A = 0, r) vs. r. The coefficients are 
only defined for r > t\ = —2.65. These results are to be compared with those in Fig. 2. 
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